Helical Fermi arcs and surface states in time-reversal invariant Weyl semimetals 
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Weyl semimetals are gapless three-dimensional topological materials where two bands touch at 
even number of points in the Brillouin zone. In this work we study a zincblende lattice model 
realizing a time-reversal invariant Weyl semimetal. The bulk dynamics is decribed by twelve helical 
Weyl nodes. Surface states form a peculiar quasi two-dimensional helical metal fundamentally 
different from the Dirac form typical for topological insulators. Allowed direction of velocity and 
spin of low-energy surface excitations are locked to the cubic symmetry axes. The studied system 
illustrates general properties of surface states in systems with common crystal symmetries. 
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Introduction- Study of gapped topological matter, 
such as quantum Hall states, topological insulators (TI) 
and topological superconductors, has become a central 
topic in condensed matter physics [lj. Topological ma- 
terials exhibit a wide range of exotic phenomena such as 
fractionalized excitations, quantized responses and un- 
usual surface states that have become signature proper- 
ties of these systems. A rapidly growing trend in search 
of new topological materials is the rising interest in gap- 
less phases. 

Weyl semimetals are three-dimensional (3d) materials 
where two energy bands touch at even number of points 
in the reciprocal space |2J. Band touching points, called 
Weyl nodes, are possible in systems where either inver- 
sion symmetry (IS) or time-reversal symmetry (TRS) is 
broken. Robustness of Weyl semimetals follows from 
topological properties of Weyl nodes acting as sources of 
the Berry curvature and defining monopoles in momen- 
tum space [2, 3J. The characteristic physical properties 
of Weyl semimetals include unusual electromagnetic re- 
sponse [4HTT] such as anisotropic DC conductivity and 
anomalous Hall effect and existence of Fermi arcs and 
topological surface states [12J and exotic broken symme- 
try states p~3| [T4] . There exists a number of promising 
theoretical proposals [12j [T5H2T] but no experimental 
verification has been achieved yet. While Weyl semimet- 
als with broken TRS have been studied extensively, the 
time-reversal invariant case has received much less atten- 
tion [22J. Weyl semimetals with unbroken TRS are not 
merely of academic interest since quite generally there ex- 
ists a Weyl semimetal phase between a TI and a trivial 
insulator phases when IS is broken [23, 24J. 

In this work we study a time-reversal invariant tight- 
binding model in the zincblende lattice with a nearest 
neighbor hopping and a spin-orbit coupled next-nearest 
neighbor hopping. The orbitals on the different face- 
centered cubic (FCC) sublattices have different onsite en- 
ergies, thus leading to IS breaking. The studied model 
is a variant of the Fu-Kane-Mele model realizing 3d TI 
phases [25J. It is remarkable that breaking IS by on- 
site potential leads to rich physics fundamentally differ- 



ent from the gapped TI phases. 

The studied model realizes a time-reversal invariant 
Weyl semimetal with twelve inequivalent Weyl nodes. 
In contrast to most of the previous studies where Weyl 
semimetals are realized by stacked 2d layers, the surface 
spectrum is identical for surfaces perpendicular to the 
three cubic axes. Populated surface states form a num- 
ber of helical time-reversed pockets bounded by Fermi 
arcs that are aligned with the cubic axes. Low-energy 
excitations in the vicinity of the Fermi arcs are highly un- 
usual consisting of particles moving only along the cubic 
symmetry axes with their spin pointing to the perpen- 
dicular direction. Thus the surface spectrum provides a 
new type of helical metal very different from typical TI 
surface states [lj and magnetic or irradiated [26J Rashba 
systems. Considering that much of the interest in time- 
reversal invariant TIs is concentrated on their helical sur- 
face metal, helical surface states on Weyl metals offer in- 
triguing alternative realization with distinct properties. 
The studied model illustrates general features of Fermi 
arc structure of time-reversal invariant Weyl semimetal 
with common point-group symmetries. 

Studied model- We are studying a variant of the tight- 
binding model on the zincblende lattice introduced by 
Fu, Kane and Mele in their pioneering work of topolog- 
ical insulators in three dimensions [25J. The zincblende 
lattice consists of two FCC lattices displaced along the 
space diagonal by one quarter of its length. The Hamil- 
tonian of the model is 

(i,j) i 

({id)) 

The first term corresponds to the nearest-neighbour hop- 
ping connecting points on different sublattices, the sec- 
ond second term represents the staggered onsite potential 
and the last term is the spin-orbit hopping between the 
second neighbours residing on the same sublattice char- 
acterized by A. We assume that the lattice sites are dis- 
tinguished by their onsite potential so that E{ = ±£ de- 



pending on which sublattice the site resides. This breaks 
IS of the diamond structure and is crucial for the appear- 
ance of the Weyl semimetal phase. The spin-orbit field is 
determined by the two bond vectors connecting the sec- 



ond neighbours as e^ 



and s = (s x . 



the Pauli matrices acting on spin. The staggered onsite 
potential breaks IS which is crucial in obtaining the Weyl 
semimetal phase. The model (fTl) gives rise to the Bloch 
Hamiltonian 



H = di(k)a x + d 2 (k)cTy + ea z 
cr z [D x (k)s x + Dy(k)s y + D z (k)s z 



(2) 



with 



d\ (fc) = t(l + cos k • ai + cos k • a 2 + cos k • a 3 ) 
d 2 (k) = t(sin k • ai + sin k • a 2 + sin k • a 3 ) 

D x (k) = A [(sink • a 2 — sink • a 3 — 
sin (k • a 2 — k • ai) + sin (k • a 3 — k • ai)]) (3) 

and (<J x ,<j y ,<j z ) are the Pauli matrices acting on sublat- 
tice space. The primitive and reciprocal vectors of the 
FCC lattice are given by 



ai 

a 2 

a 3 



(0,1,1) 
(1,0,1) 
(1,1,0) 



bi 

b 2 
b 3 



2tt 

a 

2tt 

a 

2tt 



(-1,1,1) 
(1,-1,1) 
(1,1,-1). 



(4) 



The other components D y / Z (k) are obtained by permut- 
ing lattice vectors cyclically in the expression for D x {k). 
Weyl node structure- The effective Weyl Hamiltonian 
can be found by reducing the four band model to a 
gapless and gapped two band models. This can be 
achieved by identifying the fc-dependent eigenstates of 
the spin sector ^| ± D) = ±\ ± D) where D = 
(D x (k),D y (k),D z (k)) and D = V^^A^F- In this ba- 
sis the spin part can be replaced by its eigenvalues in Q , 
leading to two band Hamiltonians for these spin sectors, 

H± = d 1 {k)cr x + d 2 (k)a y + (±D(k) + e)a z . (5) 

Assuming s > (e < 0), H+ [HJ) is always gapped 
but band touching in the H_ (-H + ) sector may occur 
if d\ — d 2 = and D = e (D = —e) for some 
k = fco- When the IS breaking potential is weaker than 
the critical value e < e c = 4|A| the band-touching con- 
ditions yield twelve inequivalent Weyl nodes in the first 
Brillouin zone at (±fc ,0,^), (0,±fe ,^), (±fe , V' )' 



,±fco,0), where sin ^ 



(0,^,±fc ),(^,0,±fc ,), (f, 

^. The Weyl node structure is illustrated in Fig. IT] 
At e = the system is gapless and the Weyl nodes of 
opposite helicity sit on top of each other. When inver- 
sion symmetry is weakly broken < |^| <C 1, the nodes 
move away from the three inequivalent X points (0, 0, — ), 




Figure 1: Red and blue circles represents the twelve 
independent Weyl nodes of opposite chirality in the first 
Brillouin zone. The node outside the zone is equivalent 

to the one in the k z = plane and illustrates how the 
nodes annihilate at critical value s = e c . 



(0,f ,0) and (f ,0,0). By increasing |^|, the Weyl 
nodes move toward the four inequivalent W points on 
the zone boundary. Since the W points on different faces 
are connected by reciprocal lattice vectors b^ they anni- 



hilate each other at the critical value 



4A I 



1, beyond 



which the system becomes a trivial insulator. The dashed 
lines illustrate this process for a pair of nodes. One way 
to understand the topological properties of the system is 
to map the Hamiltonians ([5| to the Bloch sphere. In the 
gapped phase when \e\ > e c it is not possible to cover 
the whole sphere since (±D(k) + e) has always the same 
sign when k takes all values in the Brillouin zone. Par- 
ticularly, if k is restricted to a 2d submanifold defining a 
2d Brillouin zone, the Chern number of the mapping is 
always zero for both spin sectors. The insulator phase is 
trivial without protected surface states. 

On general grounds it is known that six of the twelve 
nonequivalent Weyl nodes have positive (negative) helic- 
ity and that nodes related by TRS have the same he- 
licity. A two-fold rotational symmetry along cubic axes 
imply that the nodes connected by tt rotation have the 
same helicity. The zincblende structure also respects 
| rotations combined with a reflection with respect to 
the plane orthogonal to the rotation axis. The nodes 
related by this symmetry have opposite helicities due 
to reflection. These expectations can be verified in a 
systematic approach. In the vicinity of a Weyl node 
the effective Hamiltonian can be brought to the form 
H = v x k x a x + Vykydy + v z k z a z describing a hedgehog 
in the reciprocal space |2J. The topological invariant as- 
sociated with a Weyl point can be calculated as a flux 
through the surface containing the hedgehog and takes 
values v = sign(v x v y v z ). For a general Hamiltonian 
H = J2i ni{k)<Ji with a Weyl node at k = fco the invari- 
ant can be evaluated as v = sign [det (c^n^fc))].^. 



=fco" 



Employing this formula we can now identify the he- 
licities of Weyl points yielding ^f±k 0l o^) — —sign 6, 

V(p t ±k ,¥) = sign£ ' ^ ,±fco,o) = -signe, ^, ,±/e ) = 
sign e, u {±k0i 2^ iQ) = signs, V( ^,± ko ,) = -signs. As 
illustrated in Fig. [T] this helicity structure confirms ex- 
pectations from symmetry considerations. The gapless 
sector of Eq. ([5| describes a 3d helical semimetal with 
opposite spin textures near nodes connected by TRS. 

Fermi arcs and surface states- In the presence of a 
surface, Weyl semimetals exhibit topologically protected 
surface states. Projections of Weyl nodes to the surface 
determine end points of Fermi arcs, a line in the sur- 
face Brillouin zone separating occupied and empty states. 
The Fermi arc structure determines low-energy proper- 
ties of surface excitations. It is expected that Fermi arcs 
are formed between the nodes that annihilate each other 
at W points since a trivial insulator does not support 
surface states. 

Assuming the IS breaking parameter satisfies e > 0, 
the Weyl nodes are contained in the H_ sector with spin 
wavefunction \ — D). First we study the surface spec- 
trum when the system is terminated by a (001) plane. 
The in-plane momenta k x , k y remain good quantum num- 
bers and parametrize the spectrum in the surface Bril- 
louin zone | k x ± k y | < — . The Fermi arcs in this direc- 
tion follow from the horizontally moving Weyl nodes, the 
arcs from the vertically moving nodes shrink to points 
in the surface Brillouin zone boundary. As illustrated in 
Fig. [T] the horizontally moving nodes have lattice equiv- 
alent counterparts in the k z = ^ plane so we can con- 
centrate on the neighborhood of H_(k z = — ). In or- 
der to solve the surface modes, we define a new variable 
q = (k x ,k y , 2 ^ L ) and linearize H- in q z = k z — ^ around 
q z = which yields 

H{q) = difaK + d 2 {q)a y + [-D(q) + e] cr z + 

+ q z [d[(q)a x + d' 2 {q)a y ] + 0[q% (6) 



where d[(k) = d kz d 1 (k x ,k y ,k z ) ] 



and d' 2 {k) = 



d kz d 2 {k x ,k y ^). Defining d'{q) = ^df + d' 2 2 and 
d(q) = y/d% + d\ and employing the orthogonality 
d[(q)di(q) + d' 2 {q)d2{q) = 0, it follows that a rotation 
in the sublattice space U = e~ l9cTz ^ 2 with sin — -^ 
transforms Eq. ([6| to 

H = q z d'(q)a x + d(q)sign(k x k y )a y + [e - D(q)] a z . (7) 

The sign function determines whether (d^d^O), 
(^1,(^2,0) and (0,0,1) defines a right or left-handed or- 
thogonal set. 

By substituting q z — » — id z in Eq. (|7|), we search for 
normalizable surface states satisfying H^/ — E{k x ,k y )^ . 
Before proceeding further we must determine the appro- 
priate boundary condition that describes the vacuum- 
Weyl semimetal boundary. The system is gapped when 
e > e c and as discussed above, this phase corresponds 
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Figure 2: (a): The surface states first appear on 
separate small patches when the system is close to the 
band insulator transition f^ = 0.95 and grow together 

when the system enter deeper into the Weyl metal 

phase. White lines connecting the nodes represent the 

Fermi arcs, (b): Same as (a) but for ^- = 0.75. (c): 

0.75. 
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Energy dispersion of the surface states for 

(d): Spin texture of the surface states. On the Fermi 

arcs spin points in one of the four allowed directions. 



to a topologically trivial band insulator. Choosing coor- 
dinates so that z = coincides with the boundary, the 
surface states can found by requiring that e = £+ > e c 
for z > (outside the system) and e = £- < e c for 
z < (inside the system). From the topological point of 
view vacuum can be regarded as a band insulator with 
a very large band gap, so the above boundary condition 
is appropriate to describe a Weyl metal terminated by a 
surface |18| . 

For positive e we find a state localized on the surface 
for each (k x , k y ) satisfying D(q) > 8- given by 



ty(k x ,k y ,z) 



Ae T ^ 



y)\-D), 



(8) 



where upper (lower) signs correspond to region outside 
(inside) the system. Here | + y) denotes the state for 
which (j y | + y) = | + y) and we have included the spin 
part of the wavefunction explicitly. The momentum- 
dependent penetration depth of wavefunctions into the 
bulk is given by £± = ± _ D ( \\ £ (which is positive under 
our assumptions) . The dispersion of the surface states is 
given by 



E\kxi ky) 



sign(k x k y )d(q) 



f\j X Ob . rVyCL , v 

4tsm sin— ^— . (9) 

4 4 v J 



As plotted in Figs. [2] (a) and (b), the region for allowed 
surface states are bounded by curves where the penetra- 
tion depth £_ diverges (the black dashed curves). The 



divergence of £_ is accompanied by the merging of the 
surface spectrum to the bulk spectrum so the surface 
states cannot be distinguished from the bulk states any- 
more. The positive (negative) energy states exist in dis- 
joint, time-reversed patches of the surface Brillouin zone. 
The Fermi arcs connecting the Weyl points and satisfy- 
ing condition E(k xi k y ) =0 coincide with the coordinate 
axes k x = and k y = separating the positive and neg- 
ative energy surface pockets as illustrated in Figs. [2] (a), 
(b) and (c). Spin of the surface states, shown in Fig. [2] 
(d), is parallel to (—D x (k)^—D y (k),0)\ kz= 2 2L and have 
only an in-plane component. The surface states are he- 
lical: the in-plane momentum determines the spin state 
uniquely. However as discussed below, surface states are 
fundamentally different from the Dirac form encountered 
in many TIs. 

Allowing also negative values of £, the surface states 
exists for (k Xj k y ) satisfying D(k x ^k yj ^ L ) > \s-\. The 
spectrum and eigenfunctions are given by those presented 
above with the exception that energy changes its sign and 
spin and sublattice states are flipped | — D) — » \ + D), 
| + y) — »■ | — y). Therefore the surface spectrum is not 
continuous at e = 0. This behaviour is similar to the dis- 
continuity in the anomalous Hall insulator described by a 
Dirac equation with sign-changing mass term. Although 
the spectrum is not continuous at zero mass, the pene- 
tration depth of the surface states diverge and smoothens 
out the transition in finite systems. The solution for the 
surface spectrum by employing the linearized (long wave- 
length) Hamiltonian ^\ is valid as long as the character- 
istic spatial variation is smooth £_ ^> a. This condi- 
tion fails at the surface Brillouin zone boundary since 
d'(q) = at k y = k x ± — . However, for large enough e 
values the surface spectrum does not extend to the Bril- 
louin zone boundary and analytical approximation pro- 
vides accurate description everywhere. Also, for weakly 
broken IS f < (t^) 2 the approach fails since in that case 
the two spin sectors ([5| become nearly degenerate and 
the long- wavelength approximation cannot resolve which 
one contains the Weyl nodes. However, the case of small 
e is of limited interest since the surface states are not 
well- localized to surface anyway. 

The above analysis identifies Fermi arcs and surface 
states for (001) surfaces. A similar analysis as presented 
above can be carried out for (010) and (100) surfaces 
and as the symmetric structure of the nodes in Fig. [T] 
suggests, the spectrum for those surfaces are essentially 
identical to the spectrum of the (001) surface. 

Discussion- As discussed above, the low energy sec- 
tor of surface states consist of excitations in the vicinity 
the two pairs of time-reversed Fermi arcs. Since the elec- 
tronic dispersion is flat along the Fermi arcs, velocity 
v = dkE(k) is perpendicular to the arc and spin near 
the arcs as illustrated in Figs. [2] (c) and (d). Interest- 
ingly, velocity and spin can have only one of the four 
allowed orientations near the arcs. Therefore the allowed 



low-energy surface excitations resemble to those of two 
perpendicular Id helical metals rather than a 2d heli- 
cal metal realized on the surfaces of TI many materials. 
Time-reversed Fermi arcs also posses a perfect nesting 
property generic for ID systems. The surface state struc- 
ture could be probed by spin and angular resolved pho- 
toemission spectrum techniques or by electron tunnelling 
spectroscopy as discussed in Ref. [27J. 

One naturally wonders which properties of the studied 
model are generic and carry to other Weyl semimetals 
with unbroken TRS. The key properties determining the 
Fermi arc structure and spin textures on Fermi arcs are 
TRS and crystal symmetries. Let us assume that the 
system is terminated by a plane which admits a perpen- 
dicular axis of n fold rotations and m fold rotations com- 
bined with a reflection. From TRS we know that if there 
exists a Fermi arc between nodes at k\ and &2 there must 
exist another arc terminated by nodes at —k\ and — ki 
with identical helicities, but opposite spin. Since the sys- 
tem is invariant in n fold rotations, there must exist n 
of time-reversed Fermi arc pairs. Spin on the arc is also 
rotated but helicities of rotated nodes remain the same. 
Improper rotations also lead to m time-reversed Fermi 
arc pairs but the nodes connected by this transformation 
have opposite helicities. The spin state remains unaf- 
fected in reflection. The Fermi arc and spin structure in 
Figs. [21(b) and (c) are in agreement with the two-fold ro- 
tational symmetry and four-fold rotation combined with 
reflection. The arcs obtained by two- fold rotations and 
time-reversal coincide, so there is only two pairs of arcs. 
It is now quite straightforward to imagine allowed Fermi 
arc structures in systems with different point-group sym- 
metries. 

Summary- We studied a lattice model of a time- 
reversal invariant Weyl semimetal and its helical surface 
Fermi arcs. The surface states are identical in all three 
surfaces perpendicular to the cubic axes of the zincblende 
lattice. The low-energy excitations at Fermi arcs consist 
of helical electrons with discrete allowed velocity and spin 
orientations determined by the cubic axes. The surface 
states can be regarded as a new type of helical metal 
fundamentally different from the Dirac form realized in 
typical topological insulator surfaces. 
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